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Numerous biochemical and structural studies have 
shown that the conformation of the estrogen re- 
ceptor a (ERa) can be influenced by ligand binding. 
In turn, the conformational state of ERa affects the 
ability of the receptor to interact with a wide variety 
of protein accessory factors. To globally investi- 
gate ligand-based cofactor recruitment activities 
of ERa, we have applied a flow cytometric multi- 
plexed binding assay to determine the simulta- 
neous binding of ERa to over 50 different peptides 
derived from both known cofactor proteins and 
random peptide phage display. Using over 400 
ERa-binding compounds, we have observed that 
the multiplexed in vitro peptide-binding profiles are 
distinct for a number of compounds and that these 



profiles can predict the effect that ERa ligands 
have on various cellular activities. These cell- 
based activities include transcriptional regulation 
at an estrogen response element, MCF-7 cell pro- 
liferation, and Ishikawa endometrial cell stimula- 
tion. The majority of the compound-induced diver- 
sity in the peptide profiling assay is provided by the 
unique phage display peptides. Importantly, some 
of these peptides show a sequence relationship 
with the corepressor motif, suggesting that pep- 
tides identified via phage display might represent 
natural binding partners of ERa. These in vitro: 
cellular correlations may in part explain tissue- 
specific activities of ERa-modulating compounds. 
(Molecular Endocrinology 18: 1064-1081, 2004) 



THE ESTROGEN RECEPTOR (ER) is a multifunc- 
tional, pharmaceutical^ relevant transcription fac- 
tor that serves to transduce the activities of estrogens 
through both protein-protein and protein-DNA interac- 
tions. A wide variety of protein binding partners for ER 
have been discovered and characterized, including 
nuclear coregulator accessory proteins (1, 2), cyto- 
plasmic kinases (3, 4), and transcription factors (5, 6). 
It is clear that ER plays an important role in a number 
of interrelated signal transduction pathways. 

Abbreviations: aa, Amino acids; AF1 , activation function 1 ; 
AF2, activation function 2; AIB1 , amplified in breast cancer 1; 
AP-1, activator protein-1; cpm, radioactive counts per 
minute; DAX, dosage-sensitive sex reversal-adrenal hypopla- 
sia congenita critical region on the X chromosome; DES, 
diethylstibestrol; ERa, estrogen receptor a; ERE, estrogen 
response element; LBD, ligand-binding domain; MFI, mean 
fluorescent intensity; NCoR, nuclear hormone receptor core- 
pressor; NF-kB, nuclear factor-*B; PCA, principal compo- 
nents analysis; PGC-1a, peroxisome proliferator-activated 
receptor y coactivator-1 a; RIP140, receptor-interacting pro- 
tein 140; SERMS, selective ER modulators; SHP t short het- 
erodimer partner; SMRT, silencing mediator of retinoic acid 
and thyroid receptor; SRC-1 , steroid receptor coactivator-1 ; 
STATS, signal transducer and activator of transcription 5; 
TIF2, transcriptional intermediary factor 2. 

Molecular Endocrinology is published monthly by The 
Endocrine Society (http://www.endo-society.org), the 
foremost professional society serving the endocrine 
community. 



Two different ER genes have been identified, yield- 
ing the two sequence-related protein products, ERa 
and ER|3 (7, 8). Although ER/3 is an important partici- 
pant in the estrogen response, several studies suggest 
that ERa is responsible for much of the physiology 
associated with estrogens and estrogen replacement 
therapy (9, 10). It is commonly accepted that estro- 
gens bind to the ligand-binding domain (LBD) of ERa 
and elicit, among many other activities, an increase in 
transcriptional activity for genes containing some form 
of an estrogen response element (ERE) within the pro- 
moter region. This transactivation activity is enhanced 
by the ligand-induced association of cofactor or co- 
activator proteins, such as steroid receptor coactiva- 
tor-1 (SRC-1) or cAMP response element binding 
protein-binding protein (CBP) (11). 

Numerous crystal structures have confirmed that 
the most C-terminal a helix, typically helix 12 known as 
the activation function or AF helix, of the receptor is 
prone to conformational flexibility (12-14). Agonist li- 
gands such as estradiol and diethylstilbestrol (DES), 
tend to stabilize the AF helix into a position that ex- 
poses a hydrophobic patch or groove that can accept 
an amphipathic helix presented by cofactor proteins. 
These short helical regions, called NR boxes, usually 
occur in multiples within the cofactor protein and typ- 
ically contain an LXXLL motif (where L is leucine and X 
is any amino acid) (15). Antagonist ligands, such as 
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tamoxifen and raloxifene, force helix 12 into a position 
where it occupies the coactivator binding site and 
consequently blocks the intermolecular association of 
the coactivator amphipathic helix (12). 

The conformational flexibility of helix 12 is an impor- 
tant feature of nuclear receptors and can allow for the 
design of ER-modulating ligands that shift the dy- 
namic equilibrium between multiple states. In addition, 
nuclear magnetic resonance studies with another 
structurally related nuclear receptor have shown that 
ligands can affect the conformational dynamics in 
other regions of the LBD in addition to the AF helix 
(16). The conclusion of these structural studies was 
that subtle modifications of the ligand produce signif- 
icant effects on the conformation of the LBD in multi- 
ple regions. In light of the emerging structural and 
functional work on ERa and its LBD, there is an on- 
going quest to find new selective ER modulators 
(SERMs) that may induce unique conformations. In 
addition to the well-known molecules raloxifene and 
tamoxifen, a new class of modulators might retain 
more of the beneficial effects of estrogen without the 
undesired side effects (17). 

Similar to the conformational complexity of the LBD, 
ERa performs a myriad of signal transduction func- 
tions once bound with ligand. The complexity of ERa 
signaling is demonstrated by the observation that 
genes can be regulated either directly, through binding 
to EREs, or indirectly, through influencing the function 
of other proteins such as nuclear factor-KB (NF-kB), 
activator protein-1 (AP-1), MAPK, or STATS (6, 18-21). 
Concerning the well-studied ability of ERa to repress 
NF-KB-driven gene expression, the precise molecular 
mechanism for this activity remains unresolved and 
controversial (20). Several models have been pro- 
posed and one of these includes the direct physical 
interaction of ER with either the p50 or p65 subunit of 
NF-kB (22, 23). In addition to interaction with tran- 
scription factors, it has been shown that activated ERa 
can alter pathways in the cytoplasm, such as kinase 
phosphorylation cascades. These effects of ER are 
generally referred to as nongenomic activities. It is 
becoming more apparent that these functions are a 
significant part of the ERa activity repertoire (24). The 
development of the next generation of SERMs will 
need to consider numerous ER activities. 

To rapidly characterize the various protein-binding 
functions of ligand-activated ERa as represented by 
peptide surrogates, we have used a method that can 
essentially reconstruct a simplified in vitro cellular en- 
vironment for the receptor. This multiplexed technique 
allows simultaneous monitoring of up to 100 ligand- 
modulated interactions. Through the use of fluores- 
cently distinct microsphere populations, individual 
peptides, for example peptides representing LXXLL 
NR boxes from coactivators or sequences identified 
via phage display, can be coupled to separate micro- 
sphere populations. Flow cytometric analysis of the 
microspheres simultaneously decodes each popula- 
tion and also detects receptor binding to respective 



peptides. Use of this fluorescent microsphere multi- 
plexed binding assay has been previously demon- 
strated for nuclear receptor interactions (25), as well as 
several other applications (26, 27). 

This binding-partner profile represents, in part, op- 
portunities that the ER has to interact with partners in 
a cellular environment. In addition, phage display iden- 
tified peptide sequences have the potential to repre- 
sent currently unrecognized cellular interaction part- 
ners. These peptide interactions represent the starting 
signal for cellular outcomes driven through the ERa, 
such as gene regulation, cytokine production, cell 
proliferation, and cell stimulation. To acquire more 
evidence that peptide binding profiles reflect ERa 
biology, we have used a large number (>400) of ER- 
binding compounds that have varying degrees and 
directions of modulating activities and we have ob- 
tained both binding profiles and cell-based activities 
for these compounds. We have found that a number of 
peptide binding profiles correlate well with cellular ac- 
tivities. For this report, we have used several cell- 
based assays to explore the ability of ligand-activated 
ERa to effect 1) transcription of a reporter gene via an 
ERE, 2) the proliferation of a breast cancer cell line, 
and 3) the stimulation of an endometrial-derived cell 
line. 

As described in this paper, a simplified in vitro bind- 
ing assay for monitoring the myriad of activities of the 
ERa can be a great aid in helping to predict and 
understand the cellular outcomes of modulating li- 
gands. Using the compounds as the link between 
these two types of assays, in vitro and cellular, we 
have been able to make distinct correlations between 
peptide binding events and cell-based assay activities 
imparted by novel ERa ligands. 



RESULTS 

Peptide Tools for Compound Profiling 

Previous work has demonstrated that recombinant 
peptide libraries, such as phage display, are a power- 
ful way of identifying peptides that selectively bind to 
nuclear receptors (28-30). Our initial efforts to identify 
selective conformation-sensing peptides used a ran- 
dom 10-residue peptide library displayed as a fusion 
to the gene III coat protein of filamentous bacterio- 
phage (31). Several of the conformation-sensing pep- 
tides selected against ERa are shown in Table 1 . Typ- 
ically, with selection vs. ERa LBD bound to an agonist 
ligand such as estradiol or DES, we readily identified 
LXXLL-containing peptides after a single round of 
sorting. However, not all peptides that bound to the 
ERa LBD-estradiol complex contained the NR box 
motif. One example is the peptide denoted K7P2 as 
shown in Table 1. K7P2 binds to ERa LBD in the 
presence of many different estradiol-like ligands with 
the same behavior as a classical LXXLL, but does not 
contain a typical leucine motif. In contrast, clones 
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selected vs. antagonist-bound ERa LBD, such as 
raloxifene, tamoxifen, idoxifene, and levermeloxifene, 
yielded only non-LXXLL peptide sequences. Several 
novel ERa LBD compound-selective peptides were 
identified including I8P2 for idoxifene, R5P2 for ralox- 
ifene, and GW5P2 for GW5638 or its 4-hydroxyl me- 
tabolite, GW7604. Although we were able to success- 
fully isolate peptides with the random 10-residue 
peptide library, it is evident from sibling DNA se- 
quences of selected peptides that the size of the 
original library limits peptide diversity. 

To isolate peptides of higher affinity and greater 
specificity to the ERa LBD-GW7604 complex, we 
made a second-generation phage display library on a 
phagemid vector for soft randomization or affinity mat- 
uration (32-34). This library was based on the se- 
quence of the GW5P2 clone. Each of the ten residues 
was randomized with a bias toward the original codon 
using 80% of the bases from the original template 
sequence and an equal mixture of the three other 
bases at 20%. Three additional randomized codons 
were added to both ends of the original sequence. By 
sorting this library against ERa LBD in the presence of 
GW7604 (or other compounds with an acrylate head 
group), we identified' peptide sequences 47P3 and 
33P4. Although this peptide library was re-random- 
ized, the selected acrylate-sensing peptide sequences 
retained the original LF(E/D)RVWLRE(LyH) motif of 
GW5P2 (Table 1). An additional leucine residue was 
selected in the three randomized codon regions at 
both the N and C terminus. These extra hydrophobic 
residues may account for the increased affinity of 
47P3 for ERa compared with the original GW5P2 pep- 
tide (demonstrated by increased binding signal; see 
Fig. 2B). 

In addition to 10 peptides identified with phage dis- 
play, 48 synthetic peptides containing the binding mo- 
tifs from 18 coactivator proteins, two coregulator pro- 
teins [i.e. DAX (dosage-sensitive sex reversal-adrenal 
hypoplasia congenita critical region on the X chromo- 
some) and SHP (short heterodimer partner)] and two 
corepressor proteins (Table 1 and Ref. 25) were used. 
These peptides were typically about 25 amino acids in 
length with a biotin at the amino terminus. 

Multiplexed-Peptide Binding Assay 

To efficiently assess the effects of 405 ERa-binding 
compounds on ERa binding to 58 peptides represent- 
ing the binding motifs of coactivator, coregulator, and 
corepressor proteins as well as compound-selective 
phage display peptide sequences, we employed a 
modified microsphere-based binding assay that has 
been described previously (25). All compounds within 




1 7-(5 estradiol diethylstiibestrol 




idoxifene GW7604 



Fig. 1. Structures of ERa Tool SERM Compounds 

this set were determined to bind ERa with a pK; > 6.0. 
We included several classic, well-characterized ERa 
ligands such as estradiol, DES, raloxifene, 4-OH ta- 
moxifen, idoxifene, levormeloxifene, lasofoxifene, and 
GW7604 (35) (Fig. 1) within the 405 compound set. 
Binding of Alexa532-labeled ERa LBD is reported in 
mean fluorescence intensity units (MFI). Under the 
conditions of the assay, the MFI is proportional to 
affinity of peptide binding (25). 

In the absence of compound, ERa LBD binds to 
many of the peptides (Fig. 2A). Compounds may mod- 
ulate the interaction of ERa with peptide by enhancing 
or inhibiting this basal interaction. Typically, to dem- 
onstrate the effect of compound on ERa peptide- 
recruitment activity, we subtract the basal value from 
the value obtained in the presence of compound. 
Therefore, if a compound enhances a binding interac- 
tion, the data point will have a value greater than zero 
and alternatively, if the interaction is inhibited it will 



Fig. 2. Peptide Profile Plots for ERa 

Alexa 532-labeled ERa LBD (1 0 nwi) was incubated with 58 fluorescently unique, peptide-coupled microsphere populations for 
1 .5 h in the presence or absence of compound (1 jlim). The fluorescence associated with each microsphere was determined by 
flow cytometry. A, Binding of ERa LBD to peptides in the absence of compound. B, Binding of ERa LBD to peptides in the 
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Phage display 
peptides 




presence of compound. Each line represents the binding profile for one of 405 compounds where basal binding (no compound) 
has been subtracted. Binding profiles for tool SERM compounds are color-coded. Other compounds tested are shown as gray 
in the background. C, ERa LBD binding in the presence of tool SERM compounds {color scheme as above) for peptides selected 
by principle component analysis (basal binding has been subtracted). Most data points represent the mean from 10 separate 
binding experiments. 
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have a value less than zero. Figure 2B highlights the 
binding profiles of ERa LBD for 58 peptides in the 
presence of estradiol, DES, 4-OH tamoxifen, raloxifene, 
GW7604, idoxifene, lasofoxifene, and levormeloxifene. 
Each line represents binding data for one compound 
across the set of peptides. Other compounds exam- 
ined are shown as gray lines in the background. Es- 
tradiol and DES enhanced the binding of ERa LBD 
to coactivator and coregulator peptides as well as 
to the pan-agonist phage-identified peptide K7P2. 
In contrast, 4-OH tamoxifen, raloxifene, GW7604, 
idoxifene, and levormeloxifene all inhibited ERa LBD 
basal binding to coactivator and coregulator peptides. 
4-OH tamoxifen enhanced the binding of ERa LBD to 
the tamoxifen-selective peptide sequence a/jSV (29), 
and GW7604 promoted ERa LBD binding to the pep- 
tides GW5P2, 47P3, and 33P4. Idoxifene, and to some 
degree levormeloxifene and lasofoxifene, enhanced 
ERa LBD binding to I8P2. 

Shape-Based Clustering of Peptide 
Binding Profiles 

Many of the peptides fall into groups that either be- 
have similarly in the presence of compound or do not 
bind ERa LBD. To identify these peptide sets and 
reduce possible statistical bias, principal components 
analysis (PCA) was conducted using binding data from 
the 405 compound set. PCA reduced the number of 
peptides in the set from 58 to 1 0 representative pep- 
tides that account for approximately 99% of the vari- 
ation in the data (see Data Reduction Methods). Figure 
2C shows the same binding data from Fig. 2B after 
peptide selection by PCA. Interestingly, we have found 
that most of the coactivator binding trends for these 
compounds on ERa LBD could be described by data 
from receptor-interacting protein 140 (RIP140) (699- 
723), RIP140 (922-946), transcriptional intermediary 
factor 2 (TIF2) (676-702), peroxisome proliferator- 
activated protein y coactivator-1 a (PGC-1a) (130- 
154), and SHP (7-31). From the data shown in Fig. 2, 
B and C, it is evident that most of the compounds 
affect coactivator peptides in the same manner. In 
other words, if a compound promotes ERa LBD bind- 
ing to one LXXLL peptide it will most likely do so for 
other LXXLL peptides (especially for peptides that 
bound in the absence of compound). Subsequent 
clustering of the profile data was conducted, using the 
10 peptide variables, as shown in Figs. 3 and 4. 

Using a cutoff of 12 clusters, depicted by the den- 
drogram in Fig. 3, Fig. 4 shows the shape-based clus- 
tering results from the 405 binding profiles. A number 
of compounds yielded profiles that were relatively 
muted in terms of peptide recruitment. These com- 
pounds fall into cluster 1 (1 78 compounds) and cluster 
2 (45 compounds). Clusters 3 (three compounds), 4 
(14 compounds), and 5 (32 compounds) can be de- 
scribed as agonist or estradiol like. Compounds in 
these clusters show significant coactivator and K7P2 
peptide recruitment with little effect on ERa LBD bind- 



ing to corepressor and most phage-derived peptides. 
Clusters 6 (two compounds), 7 (33 compounds), and 8 
(59 compounds) are each unique and contain the tool 
compounds idoxifene, 4-OH tamoxifen, and ralox- 
ifene/lasofoxifene/levormeloxifene, respectively. Al- 
though ERa LBD binding to coactivator peptides is 
inhibited for all the compounds that fall within clusters 
6-8, the unique binding attribute of each of these 
clusters comes from the phage display peptides. Clus- 
ter 9 (1 compound) shows a profile somewhat similar 
to that for the 4-OH tamoxifen cluster 7. The final three 
clusters, cluster 10 (27 compounds), 11 (one com- 
pound), and 12 (10 compounds), contain variations of 
profiles that show binding to the GW5P2 peptide. In 
addition, many compounds in cluster 12 (and some 
found in cluster 10) also yield significant recruitment of 
2AP3, which is a phage peptide identified with a 
GW7604-like ligand. Overall, this analysis shows that 
different ERa ligands have the ability to differentially 
effect peptide recruitment and that these differences 
can be used to group compounds using a hierarchical 
clustering routine. 

Correlation of Cellular Function with Peptide 
Binding Profile 

Several cell-based assays were used to further char- 
acterize a subset of the 405 ERa compounds. A tran- 
sient transfection assay was used to assess the ability 
of compounds to drive reporter gene expression via 
the consensus vitellogenin ERE promoter (36). Com- 
pounds were evaluated for their ability to stimulate 
transcriptional activity (agonist mode, Fig. 5A). Those 
that did not show agonist activity were evaluated for 
their ability to antagonize or inhibit transcriptional ac- 
tivity of a subsaturating dose of estradiol (antagonist 
mode, Fig. 5B). To incorporate the pEC 50 and plC 50 
with efficacy values from the full dose-response curve 
fits, the potency measurement reported on the x-axis 
is the product of the two (where efficacy is a percent- 
age relative to either estradiol or raloxifene standards). 
In cases where we have full dose-response curves, we 
chose to display the data in this manner because it 
incorporates both measured parameters. Trends are 
the same using only efficacy values, but the correlation 
with peptide binding data is improved when account- 
ing for EC 50 or IC 50 value. As shown in Fig. 5A, a 
subset of ligands was examined for their effects on 
binding to the SRC-1(2) peptide (blue squares) and 
also for agonist activity at an ERE. These same ligands 
were also tested for binding to peptides from amplified 
in breast cancer 1 (AIB1) (607-631) (green triangles), 
K7P2 phage derived (red circles), PGC-1a (130-154) 
(yellow triangles), DAX1 (132-156) (cyan hexagons), 
and TIF2 (676-702) (red triangles). There is a general 
trend of increased ERE potency relative to coactivator 
peptide recruitment. Similarly, those compounds most 
effective at antagonizing the activity of estradiol show 
the most negative peptide recruitment values (Fig. 5B). 
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Fig. 3. Dendrogram Showing Relationship of Peptide Binding Profiles 

Peptide binding profiles were analyzed using both principle component analysis and a hierarchical clustering routine using a 
similarity measure based on a profile's shape. The resulting dendrogram shows the relationship between clusters and the vertical 
line (at 12) indicates the level chosen for display of clusters. 



A subset of the 405 compounds was also tested for 
the ability to induce cell proliferation of MCF-7 breast 
cancer cells as determined by [ 3 H]thymidine incorpo- 
ration (37, 38). A full dose-response curve was per- 
formed for these compounds and efficacies were nor- 
malized to an estradiol standard. Typically, estradiol 
yielded about an 8-fold increase in [ 3 H]thymidine in- 
corporation compared with the no compound control. 
Some compounds decreased the basal level of prolif- 
eration. For these compounds, the fit of the negative 
dose response yielded an IC 50 value and a minimum 
signal. This minimum signal was normalized to that 
obtained for raloxifene. The x-axis of Fig. 6 shows the 
potency of the compounds expressed as the product 
pEC 50 or plC 50 (where plC 50 are negative values) and 
the normalized efficacy value; the y-axis shows pep- 
tide data for these same compounds where blue 
squares represent binding data for SRC-1 (2) (676-700) 



peptide and other colored symbols represent five ad- 
ditional peptides. As shown, there is a general trend 
toward increased coactivator recruitment and in- 
creased ability to cause MCF-7 cell proliferation. 

The uterotrophic activity of compounds was assessed 
with an Ishikawa (a human endometrial carcinoma cell 
line) cell stimulation assay using intrinsic alkaline phos- 
phatase activity as a reporter (39, 40). Maximum alkaline 
phosphatase activity values were obtained for each 
compound (typically at concentrations greater than 1 
/xm), and these values were normalized to an estradiol 
control. For example, 4-OH tamoxifen, raloxifene, and 
GW7604 typically gave %E max values of 14, 7, and 4 (41). 
Figure 7A shows Ishikawa data for a subset of com- 
pounds plotted vs. SRC-1 (2) peptide binding (blue 
squares) for the same set of compounds. Data for five 
other peptides are also shown [other colored symbols). 
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Fig. 4. Categorization of ERa-Binding Compounds by Clustering of ERa LBD Peptide Binding Profiles 

Profile plots of ERa LBD binding to peptides in the presence of various compounds are shown. Peptides shown are those 
chosen via the PCA-based method described in Data Reduction Methods. Profiles are grouped by shape-based cluster ID, which 
is represented in Fig. 3. Each line represents binding data for one compound where ERa LBD binding in the absence of compound 
has been subtracted. Each cluster contains peptide-binding data for one or more ERa compounds. The labels described in Fig. 
2 identify profiles for tool SERM compounds. 



Figure 7A shows that there is a general correlation be- 
tween coactivator peptide recruitment and ability to 
stimulate Ishikawa cells. Figure 7B shows the relation- 
ship between Ishikawa stimulation and recruitment of 
three phage peptides, a//3V (tamoxifen-selective; pink 
triangles), I8P2 (idoxifene-selective; cyan circles), and 
GW5P2 (GW7604-selective; yellow squares). As shown, 
for compounds that showed significant recruitment of 
either the odfiV or I8P2 peptides, Ishikawa %E max values 
were less than 20%. For compounds that recruit the 
GW5P2 peptide, the Ishikawa %E max values were even 
lower (<5%). 

In Figs. 5-7, we have shown that there is a general 
correlation between coactivator/phage peptide bind- 
ing and the modulation of ER-driven cell-based activ- 
ities imparted by these compounds. To further explore 
this relationship, Fig. 8 displays cell-based measure- 
ments for ERE, MCF-7, and ishikawa assays in the 
same clusters as defined by the peptide binding data 
represented in Fig. 4. Figure 8 shows that, in general, 



compounds that cluster together based on peptide 
profile also show similar cell-based activities. For sev- 
eral of the clusters, the cell-based data show highly 
consistent results within the cluster. For example, 
cluster 3 (with estradiol and DES) shows the expected 
high Ishikawa stimulation, MCF-7 proliferation and 
ERE activity. Cluster 7, containing 4-OH tamoxifen, 
shows low Ishikawa stimulation, MCF-7 proliferation 
and ERE activity. Compounds in cluster 10, which 
contains GW7604, show very low Ishikawa %E max 
values (between 2 and 5%). Clusters 1 1 and 12, which 
like cluster 10 is defined by strong binding to the 
GW5P2-like peptides, also show a trend for low Ish- 
ikawa response. Compounds in cluster 8, which in- 
cludes raloxifene, lasofoxifene, and levormeloxifene, 
yield relatively low Ishikawa %E max values (between 5 
and 20%) and show the ability to inhibit both MCF-7 
proliferation and ERE activation. 

Interestingly, for the compounds in clusters 1 and 2, 
where the peptide profiles were relatively flat, there did 
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ERE (pEG50)(%max) 




ERE (plC50)(%max) 

Fig. 5. Effects of Compounds on ERE-Driven Reporter Gene 
Transcription as compared with ER« LBD Binding to Select 
Cofactor Peptides 

ERa LBD peptide binding data were determined as de- 
scribed in Fig. 2. Each point represents binding data for one 
compound where ERa LBD binding in the absence of com- 
pound has been subtracted. Data for six peptides are repre- 
sented by the following colors: SRC-1(2) (676-700) (blue 
squares), AIB1 (607-631) (green triangles), K7P2 phage de- 
rived (red circles), PGC-1« (130-154) (yellow triangles), DAX1 
(132-156) (cyan hexagons), TIF2 (676-702) (red triangles). 
ERE data for the same compounds are plotted on the x-axis. 
The ERE assay was run in either (A) agonist or (B) antagonist 
mode as described in Materials and Methods. The % max 
values were determined by normalizing to either estradiol or 
raloxifene, respectively. The solid line indicates the best fit of 
all data points. The dotted line indicates 95% confidence 
interval and the dashed line is 95% prediction interval. 



not appear to be much relationship among the cell- 
based data. Essentially, it appears that where there are 
no extreme positive or negative peptide descriptors for 
these compounds, these ligands can have a variable 
response in cells. One explanation for this is that our 
peptide interaction set does not contain enough inter- 
actors or descriptors to clearly differentiate com- 
pounds within these clusters. Also, cluster 5 shows 
heterogeneous cell-based responses. One attribute 
that clusters 1, 2, 5, and 8 have in common is the 
absence of a phage display peptide(s) to characterize 
the group. It is possible that additional phage-derived 
conformation sensing peptides could help in subcat- 
egorizing these clusters. 




.1000; .400' 0 500 



MCF-7 Potency [{pECso or p!C50X%cpm)] 

Fig. 6. Effects of Compounds on MCF-7 Cell Proliferation 
as Compared with ERa LBD Binding to Select Cofactor 
Peptides 

ERa LBD peptide binding data were determined as de- 
scribed in Fig. 2. Each point represents binding data for one 
compound where ERa LBD binding in the absence of com- 
pound has been subtracted. Data for six peptides are repre- 
sented by the following colors: SRC-1(2) (676-700) (blue 
squares), AIB1 (607-631) (green triangles), K7P2 phage de- 
rived (red circles), PGC-1 a (130-154) (yellow triangles), DAX1 
(132-156) (cyan hexagons), TIF2 (676-702) (red triangles). 
MCF-7 cell proliferation data for the same set of compounds 
are shown as an aggregate value incorporating both the 
pEC S0 or IC 50 and efficacy measurements. For %cpm deter- 
minations, data for all compounds were normalized to either 
estradiol (max) or raloxifene (min). 



Correlation of Cell-Based Compound Effects with 
Corepressor Binding 

Several recent studies have examined the role of an 
ERaxorepressor interaction in mediating the tissue- 
selective effects of tamoxifen and raloxifene (42-44). 
We used the 405 compound set to search for ligands 
that might enhance the interaction between ERa LBD 
and a peptide derived from either nuclear hormone 
receptor corepressor (NCoR) or silencing mediator of 
retinoic acid and thyroid receptor (SMRT). As shown in 
Fig. 9, A and B, a number of compounds were shown 
to enhance the recruitment of both the NCoR 1 [amino 
acids (aa) 2040-2065] and SMRT 2 (aa 2329-2354) 
peptides. The compounds that were most successful 
at enhancing this interaction are those with an acrylate 
or carboxylate head group extension (like GW7604; 
shown as green squares). In addition, some novel 
head groups were also found to induce the same 
response and cluster with the acrylate-containing 
compounds (red squares). As shown in the two exam- 
ples in Fig. 9, the same compounds that increased 
corepressor peptide binding also increased binding to 
the GW5P2 phage display peptide and the GW5P2- 
iike peptide, 47P3. 

Examination of the peptide sequences found with 
the ERa:GW7604 complex revealed that these pep- 
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Fig. 7. Effects of Compounds on Ishikawa Cell Stimulation 
as Compared with ERa LBD Binding to Select Cof actor 
Peptides 

ERa LBD peptide binding data were determined as de- 
scribed in Fig. 2. Each set of points represents binding data 
for one compound where ERa LBD binding in the absence of 
compound has been subtracted, Ishikawa cell stimulation 
data are shown for the same set of compounds. For the 
Ishikawa %E max value, each compound was normalized to a 
standard estradiol stimulation curve (100%). A, SRC-1(2) 
(676-700) (blue squares), AIB1 (607-631) {green triangles), 
K7P2 phage derived (red circles), PGC-1« (130-154) (yellow 
tn'angles), DAX1 (132-156) (cyan hexagons), TIF2 (676-702) 
(red triangles). B, GW5P2 phage display peptide (yellow 
squares), I8P2 phage display peptide [cyan circles), oJfN 
phage display peptide (pink triangles). 



tides share the same hydrophobic residue spacing as 
corepressor motifs (Fig. 1 0). This alignment shows that 
the typical corepressor motif of LXX(H/I)IXXX(L/I) is 
represented in the majority of the acrylate-sensing 
phage display peptides. In addition, a Leu at the -1 
position and the basic amino acid at the +7 position 
were selected in both the primary random 10-residue 
peptide library and the second-generation GW5P2- 
based library. Selection from the latter library identified 
a Thr at the -2 position and a Leu at the -3 and +12 
position to further extend the presumed amphipathic 
helix. Also, Phe and Trp were preferentially selected at 
the +1 and +5 positions, where the corepressor con- 
tains the smaller hydrophobic amino acids Leu and lie, 
respectively. 



DISCUSSION 

The search for novel ERa ligands has previously relied 
heavily on cell-based predictive assays. A number of 
studies, using protease digestion (45), crystallography 
(12, 46), and phage display (28, 29), have elegantly 
shed light on the molecular mechanism for ligand- 
induced receptor conformations for ERa. Such studies 
have highlighted the importance of receptor dynamics 
in governing outcomes of ER-binding compounds. An 
important step in further understanding the details of 
ERa signal transduction and the consequent discov- 
ery of new modulating ligands is further defining the 
link between receptor conformation and biological ac- 
tivity of SERMS. Here, we have used various ERa- 
binding small molecules as tools to probe both re- 
ceptor conformation, via peptide interactions, and 
cell-based activity. In general, these results show that 
peptide interaction profiles correlate well with, and 
could be used to predict, the biological consequence 
of ligand binding within a cell as described by a set of 
cell-based predictive assays. These correlations fur- 
ther highlight the dominant role that ligand-induced 
conformation and association with cofactor have on 
ERa cellular function. 

There has been tremendous effort toward under- 
standing the differential effects of SERMs and the 
effects that they cause at the cellular level (1 7). One 
approach has been to use differential gene expression 
to fingerprint and classify ER ligands based on their 
ability to alter RNA transcript levels (47, 48). It is be- 
coming clear from these types of studies that gene 
expression patterns help explain why some ER-mod- 
ulating ligands display both beneficial outcomes and 
undesired activities such as activity in uterine or breast 
tissue. Future work should allow assigning specific 
genes to particular biological effects elicited by 
SERMs, such as retention of bone mineral density (49). 
However, parsing through the enormous numbers of 
available and newly synthesized ER-binding com- 
pounds at the cellular level is a daunting task. The 
results presented here demonstrate that the ligand- 
induced conformation of ER can be probed through 
protein-peptide interactions and these binding profiles 
may be used to characterize and bin new SERMs to 
aid selection of novel compounds. Compounds that 
yield profiles unlike those of known SERMs can be 
rapidly identified and flagged for further cell-based or 
animal-based work. This relatively high-throughput 
technology should help guide ER and other nuclear 
receptor chemistries for fine-tuning modulators of the 
future. 

Interestingly, we have made these in vitro to cellular 
correlations even though our multiplexed peptide 
binding assay uses purified LBD. In fact, our ERa 
expression construct contains an abbreviated form of 
the LBD (aa 299-554) that deletes the extended strand 
after the AF helix (F domain). Although we observe 
excellent correlations, expanding efforts to include 
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Fig. 8. Relationship of Cellular Data to ERa LBD Peptide Binding Profiles for ERa-Binding Compounds by Clustering 

Cell-based assay data for Ishikawa, MCF-7, and ERE are shown in each cluster as organized by peptide profile characteristics. 
Cell-based data were generated as described in Figs. 5-7 and Material and Methods. Clusters 1-12 were determined based on 
the ER or LBD peptide profiles described in Fig. 4A. Each point represents eel I- based data for one compound and each cluster 
is based on peptide-binding data for one or more compounds. Cell-based data for tool SERM compounds are coded as follows: 
estradiol (blue-green), DES {dark blue), 4-OH tamoxifen (red), GW7604 (black), idoxifene (green), levormeloxifene (light blue), 
lasofoxifene (blue), and raloxifene (pink). Other compounds tested are shown as gray open circles in the background. The scale 



for the celi-based values is as follows: Ishtkawa E% n 
0 - 1000; ERE pEC 50 x % max , 0 - 1000; ERE plC 5) 



1X , 0 - 100%; MCF-7 pEC 50 x %cpm, 0 - 1000; MCF-7 plC 50 x %cpm, 
x % ma „ 0 - 1000. 



full-length receptor and peptide/protein/DNA binding 
partners that recognize other regions of ER (such as 
the AF-1/DNA-binding domain/hinge regions) should 
increase the amount of information garnered. It is well 
documented that ERa has numerous activities driven 
through the N-terminal AF1 domain. For example, it 
has been shown that the AF1 domain is phosphory- 
lated upon ligand binding and that other ligand- 
induced activities can occur via cytoplasmic signaling 
events (50). The relationships that we derive in this 
current study demonstrate that ligand-influenced AF2 
function is a dominant feature of ERa signal transduc- 
tion, and it is likely that AF1 phosphorylation and its 
activities are linked to the AF2 in the context of the 
full-length receptor. Although our data suggest that 
most of the receptor function correlates with ligand- 
dependent AF2 function, it is important to consider the 
possibility that in vitro profiling may not completely 
correlate with ERa LBD binding and function. A case 



could be made that nongenomic effects, possibly 
driven through an alternate signal transduction path- 
way such as recently shown for progestin (51), may 
yield discrepancies between LBD profiling and cell- 
based activities for ERa ligands. Additionally, lack of 
correlation between peptide binding characteristics in 
vitro and the effect of compounds on a cell could be 
due to lack of cell penetration by a compound or due 
to metabolic alteration of a compound upon entering a 
cellular environment. 

Although a number of studies have suggested dif- 
ferential coactivator recruitment by ERa ligands (52, 
53), we have found, under our experimental condi- 
tions, that binding characteristics for most cofactor- 
derived LXXLL peptides track with one another. Gen- 
erally, all estradiol-like compounds proportionally 
increased LXXLL peptide binding relative to basal 
binding and all tamoxifen- or raloxifene-like com- 
pounds proportionally decreased LXXLL binding. Ap- 
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Fig. 9. Comparison of Compound Effects on ERa LBD Binding to Phage Display and Corepressor Peptides 

Binding of compound-bound ERa to the respective peptide was measured using fluorescent microspheres and flow cytometry 
as described in Materials and Methods and Fig. 2. Each point represents binding data for one compound where ERa LBD binding 
in the absence of compound has been subtracted. The three colors represent different classes of compounds: green depicts 
compounds that contain either an acrylate or carboxylate head group (GW7604 is one example), biue depicts compounds that 
contain a tertiary amine head group (4-OH tamoxifen is one example), and red shows those compounds with a unique head group. 
A variety of different core scaffolds are represented in addition to the triphenylethylene. 



plication of PCA deselected 39 of the 44 cofactor 
peptides used in this study. In other words, the five 
cofactor peptides shown in Fig. 2C were adequate to 
fully describe the compound-induced binding trends 
of all cofactor peptides. It is possible that very subtle 
differences in cofactor recruitment with various com- 



pounds may require full-length proteins, where sites 
distal to the coactivator cleft (on the receptor) and/or 
NR box (on the cofactor) may play a role. Clearly 
though, as noted by previous studies (54, 55), there 
are distinct LXXLL preferences for ERa LBD where, as 
shown in Fig. 2, A and B, peptides from RIP140, 
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Fig. 10. Alignment of Phage Display Peptide Sequences with NCoR and SMRT Interaction Motifs 

The peptides identified using random peptide phage display against ERa complexed with GW7604 or other acrylate head group 
compounds are shown aligned with the corepressor binding motifs from both NCoR and SMRT. The phage display peptides are 
GW5P2, 2AP3, 47P3, 33P4, and 7)3-16 [which was identified by Connor ef al. (59)]. 



SRC-1, TIF2, AIB1, PGC-1a, SHP, and DAX yield the 
highest binding signal. It is interesting to note that four 
of these five peptides, in addition to the phage-derived 
peptide K7P2 (LLRWVL), show a weak sequence re- 
lationship of <j> LKXLL, where <f> represents a hydro- 
phobic side-chain (Leu, Val, or He) and K represents 
Lys in addition to the other basic residues Arg, Gin, 
and His (Table 1) (56). The DAX peptides deviate 
slightly by containing a Tyr residue at the +2 position 
(ILYXM/LL). These peptide sequences, notably K7P2, 
DAX (66-90), and DAX (132-156), also suggest that 
ERa is able to tolerate either a Met or Tyr in the +4 
position. A full understanding of the detailed molecular 
recognition aspects of the preferential binding of cer- 
tain peptides will require further structural and func- 
tional analysis. 

Figures 5-7 show that there is a consistent relation- 
ship between compound-induced coactivator peptide 
recruitment and the effect of compound on ERa- 
driven cell response. To more accurately show this 
correlation, we have assembled various cell-based ac- 
tivities for a subset of the tool compounds and orga- 
nized these data into clusters based on similarity of 
peptide profiles (Fig. 8). As shown, most of the cell- 
based data are similar in each particular cluster. The 
most obvious exceptions are data shown in cluster 1 . 
The peptide profiles for this cluster are relatively 
unchanged from the peptide profile of ERa in the 
absence of compound (i.e. the change in relative 
fluorescent intensity is approximately zero for all pep- 
tides). Although the compounds in this cluster bind 
with relatively high affinity, they do not appear to in- 
duce a unique, stabilized conformation. These types of 
ligands may act as competitive antagonists by keep- 
ing the receptor in an apo-like state. In this situation, 
the ligand-bound LBD structure would be predomi- 
nantly a dynamic state similar to the ligandless recep- 
tor. It will be interesting to pursue whether or not 
peptides or other binding partners will be able to fur- 
ther discriminate compounds that appear to induce a 



somewhat neutral conformation. With additional pep- 
tides, the correlation of peptide profile with cellular 
activity for ligands within this cluster may improve. 

It has been noted that tamoxifen, although devel- 
oped as an antiestrogen for the treatment of breast 
cancer, has estrogenic activities in certain tissues (57). 
In particular, tamoxifen exhibits minimal, but signifi- 
cant uterotrophic activity relative to raloxifene. With 
respect to the LXXLL cofactor peptides, both tamox- 
ifen and raloxifene induce very similar peptide profiles 
(Fig. 2B). However, the binding profiles are clearly 
distinct with respect to the tamoxifen-selective pep- 
tide a/j8V (29) (Fig. 2C). Interestingly, compounds that 
induce recruitment of the peptide also yield an Ish- 
ikawa stimulation value less than 20% (Fig. 7B). One 
explanation for this relationship may be that the 
tamoxifen-selective peptide represents an uncharac- 
terized cofactor that may be expressed in uterine 
tissue but not breast tissue. Although it has been 
demonstrated by chromatin immunoprecipitation ex- 
periments (44), we were not able to detect a significant 
difference with SRC-1 peptide recruitment between 
tamoxifen and raloxifene in vitro. This may be related 
to a lack of an AF1 region on ERa or lack of full-length 
SRC-1 cofactor as regions outside the LXXLL may 
play a partial role in recruitment. In addition, we were 
not able to detect a significant interaction between 
ERa and corepressor peptide with either tamoxifen or 
raloxifene. 

To even further discriminate compounds with low 
uterine stimulatory activity via peptide profiling, we 
have selected peptides specific for two compounds, 
idoxifene and GW5638, or its 4-hydroxyl version 
GW7604. These compounds yield Ishikawa stimula- 
tion values less than 15% relative to estradiol. As 
shown in Fig. 2, B and C, the two selected peptides, 
I8P2 and GW5P2, are specifically recruited to ERa 
upon binding of idoxifene or GW5638/GW7604, re- 
spectively. In addition, our test compound set con- 
tains several molecules similar to idoxifene and 
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GW5638/GW7604 and these compounds also induce 
recruitment of the selective peptides (Fig. 2C). Similar 
to the a//3V tamoxifen-selective peptide, comparison 
of I8P2 and GW5P2 peptide recruitment with Ishikawa 
cell stimulation, shows that compounds that induce 
ERa binding to these peptides also exhibit uterine cell 
activity approximately less than 20% (Fig. 7B). From 
the comparison of peptide and Ishikawa activities, the 
GW5P2 peptide appears to predict an even lower uter- 
ine activity. Although on!y several compounds are 
shown, we have examined more than 15 compounds 
that induce recruitment of the GW5P2 peptide and all 
these compounds display an Ishikawa activity less 
than 5% relative to estradiol (data not shown). 

Sequence alignment shows that there is a relation- 
ship between the set of phage peptides isolated 
against the ERa:GW5638/GW7604 complex and the 
corepressor peptides, from both NCoR and SMRT 
(Fig. 10). Compound-induced recruitment shows that 
both these phage peptides and the corepressor pep- 
tides are recruited by same compounds in a linear 
relationship (Fig. 9, A and B). This suggests that these 
peptides may interact at the same site on ERa and 
have the same mode of interaction with ERa. It is 
possible that the GW5P2 and 47P3 peptides identified 
by phage display are affinity-optimized versions of 
corepressor peptides. If this is the case, the hydro- 
phobic face of the amphipathic helix would be approx- 
imately one helical turn longer than the typical coac- 
tivator peptide and the interaction on the receptor 
would require repositioning of the AF helix as shown in 
the recent peroxisome proliferator-activated recep- 
tora/SMRT peptide crystal structure (58). As shown in 
Fig. 9, the binding signal for the GW5P2 and 47P3 
phage peptides is approximately 3- to 5-fold higher 
than that for the corepressor peptides. One possibility 
for the increased signal (due to increased affinity) is 
that the selected Phe at +1 , Val at +4, and Trp at +5 
provide more hydrophobic surface area for packing 
against the extended corepressor cleft. It is interesting 
that the acrylate head group compounds optimally 
create this molecular recognition surface on ERa. Us- 
ing an LXXLL focused library, Connor et al. (59) also 
identified a similar peptide consensus from phage dis- 
play results with ERa:GW7604. Further structural and 
functional work will be required to elucidate the pre- 
cise molecular recognition of these peptides and rel- 
evance of their similarity to the corepressor motif. 
From this example, and the fact that phage display 
peptides contribute a large degree of the prediction 
information for peptide profiling, it appears that the 
selected peptides might represent naturally occurring 
protein partners, some of which may remain to be 
discovered. 

As mentioned briefly above, the compound-selec- 
tive peptides offer a great deal of diversity to the 
peptide-binding profiles in this type of approach. For 
examining new ligands, the profiles can give clues that 
a new ligand is unique compared with tool com- 
pounds, even if a specific peptide for that new ligand 



is not in the peptide set. In fact, compounds that 
induce either novel or neutral profiles can be selected 
for further random peptide work to identify new pep- 
tides to add to the screening set. In contrast to a static 
picture provided by a crystal structure, this type of 
profiling analysis provides a sense of the activities of 
the receptor, which is a balance between various "on" 
and "off" states. In effect, the peptide interaction pro- 
file is sensitive to the dynamics of the receptor when 
bound with ligand. 

Much of the detail for how ligand-induced receptor 
dynamics translates to signal transduction capacity in 
a cellular environment has yet to be completely de- 
fined. Our data add to the mounting evidence that the 
allosteric coregulator recruitment activity of ERa, as 
well as other members of the nuclear receptor family, 
is a dominant feature of receptor function that ulti- 
mately links biological activity to the ligand. Visualiza- 
tion of the multitude of ligand-induced peptide profiles 
further emphasizes the structural and functional plas- 
ticity of the nuclear receptor LBD. Understanding how 
to manipulate the numerous pathway activities of ERa 
with ligands, even with many current unknown factors, 
remains a goal for both drug and tool compound dis- 
covery. This microsphere-based compound profiling 
strategy, coupled with random peptide library technol- 
ogies, can be a valuable tool to rapidly assess com- 
pound-induced conformations of ERa and can aid in 
the early identification of unique potential drug 
candidates. 



MATERIALS AND METHODS 



Protein 

The ligand binding domain of ERa was expressed in Esche- 
richia co//, purified to homogeneity, and biotinylated similar to 
that described previously for ERj3 LBD (25). The expression 
construct consisted of aa 299-554 and mass spectrometry 
was used to verify the mass of the purified protein. Full-length 
ERa was purchased from Invitrogen (Carlsbad, CA). 

Selection of ERa-Binding Compounds 

All compounds were either synthesized in the GlaxoSmith- 
Kline Medicinal Chemistry Laboratories or were purchased 
from commercial sources. Solid compounds were dissolved 
in dimethylsulfoxide (10 itim) and subsequently diluted with 
buffer before use. An in vitro [ 3 H]estradiol displacement assay 
was used to select ERa-binding compounds from a large set 
of candidate compounds (data not shown). The compounds 
selected for further analysis in the microsphere-based pep- 
tide binding assay and the cell-based assays had pK/s > 6. 
The average pKj for the 405 compounds in this test set was 
7.2 with sd of 0.8. The compound set included several clas- 
sic, well-characterized ERa ligands such as estradiol, DES, 
raloxifene, 4-OH tamoxifen, idoxifene, levormeloxtfene, laso- 
foxifene, and GW7604 (35) (Fig. 1). 

Phage Display 

A modified polyvalent phage display vector, fGWg3, was 
constructed using the original vector fTC (60). fGWg3 con- 
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tains a flag tag epitope (AspTyrLysAspAspAspLys) in frame 
immediately upstream of the sequence for the Ml 3 genelll 
protein. This vector was also engineered to contain a dual 
Bbs\ site between the signal sequence of the M13 gene ill 
and the flag tag. A random peptide library, fGWX10, was 
constructed by inserting a DNA cassette containing ten con- 
secutive NNK codons, where N represents an equal molar 
ratio of A, C, G, and T and K represents G and T, into Bbs\ 
digested fGWg3. This library displays a random 10-residue 
peptide sequence with flanking regions as follows: NH2- 
EDGGSX 10 (GGGGS) 3 -g!ll protein. The DNA cassette encod- 
ing the random peptide sequence was constructed using a 
synthetic ligase chain reaction procedure (61) with the fol- 
lowing four oligonucleotides: 5'-TGGGTCCIII flllllllll llllillill 
IIIIIIIGGA GGCGGGGGTA-3', 5'-GCGGCGGTGG AGGGT- 
CTGGG GGAGGTGGCT CTGGGAAGAC CT-3\ 5'-TTCC- 
AGGTCT TCCCAGAGCC ACCTCCCCCA GACCC-3', and 
5'-TCCACCGCCG CTACCCCCGC CTCCMNNMNN MNNM- 
NNMNNM NNMNNMNNMN NMNNGGA-3', where I repre- 
sents inosine, M represents equal molar of C and A, N rep- 
resents equal molar of A, C, G, and T. E. coli MC1061 cell 
electroporation and phage library preparation have been de- 
scribed previously (62). After ligation and transformation, a 
small fraction of the library was plated on LB/tetracycline 
plates. The library contained approximately 10 10 transfor- 



mants and individual clones were analyzed by both DNA 
sequencing and phage ELISA using immobilized anti-Flag 
monoclonal antibody. 

A second-generation phagemid display library was con- 
structed based on the selected peptide sequence for ERa 
plus GW5638 (Table 1). The library oligonucleotide was made 
to re-randomize the nucleotide sequence encoding the 
GW5638 selected peptide. An oligonucleotide for cassette 
mutagenesis was synthesized with base mixtures where each 
base position was fixed at 80% based on the original tem- 
plate sequence and an equal mixture of the three other bases 
was used at 20%. In addition, three NNS, where S represents 
an equal molar ratio of G and C, randomized codons were 
added to both ends of the re-randomized oligo segment. This 
oligonucleotide mixture sequence served as a template for 
PCR amplification using the following sequences as primers: 
5'-ccg gcc atg gcc gaa gac ggt ggg-3' and 5'-tac cag gtt cca 
ggt ctt ccc aga gcc acc tec ccc aga ccc tec acc gcc get acc 
ccc gcc tcc-3'. The PCR amplified cassette was ligated into 
a Bbsl-digested, gel purified pGW BBS6g3 phagemid vector. 
The ligated library was electroporated into XL1 -Blue cells and 
titered to yield 1.25 x 10 9 individual transformants. Random 
clones were selected and analyzed by DNA sequencing to 
show that 94% of the clones had correct inserts containing 
randomized DNA sequences. 



Table 1 . Peptide Sequences Used to Profile ERa-Modulating Compounds 



Peptide 


Amino Acid Residues* 


Sequence 


Phage Display Peptides 






K7P2 


NA 


EDGGSLLRWYLEHEFGGGGSGK (Ahx-B) h - C0NH 2 


o/jSV* 


NA 


SSPGSREWFKDMLSRGGGGSK (Ahx-B) -CONH 2 


I8P2 


NA 


EDGGSWVEWVEEERRGGGGK (Ahx-B) -CONH 2 


R5P2 


NA 


EDGGSRRFTFQYGSV (GGGGS) 3 GK (Ahx-B) -C0NH 2 


2AP3 


NA 


EDGGSPLWQIFSRELGGGGSK (Ahx-B) -C0NH 2 


GW5P2 


NA 


EDGGSLFERVWLREL (GGGGS) 3 GK (Ahx-B) -CONH 2 


47P3 


NA 


QPAMAEDGGSALTLFDRVWLREHQPLGGGGSK (Ahx-B) -CONH 2 


33P4 


NA 


QPAMAEDGGSVLTLFERVWLRELSCLGGGGSK (Ahx-B) -C0NH 2 


70-1 6 d 


NA 


ASRHFLINQHLYKLLQDTDIWSRLGGGGSK (Ahx-B) -C0NH 2 


Coactivator Peptides® 






RIP140 


699-723 


B - SGSEI ENLLERRTVLQLLLGNPTKG - CONH 2 


RIP140 


922-946 


B - EHRS WARES KS FNVL KQLLL S ENC V - C0NH 2 


TIF2 


676-702 


B -GSTHGTSLKEKHKILHRLLQDSSSPVD- CONH 2 


TRAP 220 


590-614 


B-GHGEDFSKVSQNPILTSLLQITGNG-CONH 2 


TRAP 220 


631-657 


B - PVS SMAGNTKNHPMLMNLLKDNPAQ - CONH 2 


RAP250, PRIP 


874^898 


B - F PVNKDVTLTS PLL VNLLQSD I SAG - C0NH 2 


RAP250, PRIP 


1476-1500 


B-NKNLVSPAMREAPTSLSQLLDNSGA-CONH 2 


PPAR 7 CoA-1 a 


130-154 


B-DGTPPPQEAEEPSLLKKLLLAPANT- C0NH 2 


Coregulator Peptides 






SHP 


7-31 


B - AC PCQGAAS RPA I L YALL S S S LKAV - C0NH 2 


SHP 


104-128 


B-VTFEVAEAPVPSILKKILLEEPSSS-CONH 2 


DAX 


1-23 


B - MAG ENHQWQG S I L YNMLMS AKQT - CONH 2 


DAX 


66-90 


B - CCFCGKDHPRQGS ILYSMLTSAKQT - C0NH 2 


DAX 


132-156 


B - CCFCGEDHPRQGSILYSLLTSSKQT- CONH 2 


Corepressor Peptides 






NCoR 1 


2040-2065 


B - GQ V P RTHRL I TL ADH I C Q I I TQD FAR - CONH 2 


NCoR2 


2251-2275 


B-GHSFADPASNLGLEDIIRKALMGSF-C0NH 2 


SMRT1 


2124-2149 


B - PGVKGHQRWTLAQHI S EVITQDYTR - CONH 2 


SMRT2 


2329-2354 


B - S Q A VQ EHAS TNMGL E AI IRKALMGKY - C0NH 2 



a As specified in GenBank. 

h Ahx-B represents 6-aminohexanoic acid-biotin conjugate. 

c Sequence reported in Ref. 29. 

d Sequence reported in Ref. 59. 

a See Ref. 25 for other coactivator peptide sequences. 

NA, Not applicable. Boldfaced letters in the peptide sequences represent consensus residues in coactivator and corepressor 
peptides. 
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In preparation for selection of specific peptide clones, 500 
nM of either biotinylated ERa LBD or biotinylated ERa FL in 
the presence of 1 ixm compound of interest was immobilized 
on a neutravidin-coated microtiter plate well. Protein was 
immobilized for at least 1 h at room temperature and then 
wells were washed three times with PBS. For the first round 
of selection, approximately 10 12 phage particles were added 
in casein blocking buffer (0.89 g dibasic sodium phosphate, 
2.5 g Hammarsten casein, and 2.9 g sodium chloride in 250 
ml ddH 2 0) plus the compound of interest. Phage were al- 
lowed to bind for 2 h at room temperature. After washing with 
PBS/0.025% Tween 20 to remove unbound phage, bound 
phage were eluted with 0.1 n HCI plus 1 mg/ml BSA, the 
solution in the well was adjusted to pH 2.2 with glycine and 
then neutralized with 2 m Tris-HCI. An aliquot of the eluted 
and neutralized phage was added to 1 ml of confluent K91 
cells. Cells were infected with phage particles for 1 5 min at 37 
C, and then cultures were grown overnight at 37 C. The 
following day, phage were harvested by standard methods 
using 10% polyethylene glycol-8000/2.5 m NaCI in prepara- 
tion for the next round of selection. 

After each round of phage sorting, eluted phage were 
titered to determine population enrichment. Typically after 
two to three rounds of selection, enrichment was observed 
and individual clones were analyzed by DNA sequencing. 
Phage ELISA (63) was used to assess both peptide affinity for 
immobilized ERa and the selectivity of the peptide for the 
receptor in the presence of ligand. Based on phage ELISA 
signal and selectivity vs. a panel of known SERMs, individual 
peptides were chosen for synthesis. 

Synthetic Peptides 

Biotin-containing peptide sequences derived from coactiva- 
tor, corepressor, and coregulator proteins or identified by 
phage display were synthesized by the Synpep Corp. (Dublin, 
CA), or American Peptide Co. (Sunnyvale, CA), Each peptide 
was purified by reverse- phase HPLC. Amino acid composi- 
tion and peptide concentration were determined by conven- 
tional methods for each peptide sequence. Typically, these 
peptides were biotinylated during synthesis at the e-amino 
group of lysine of the G 4 SK C terminus. This orientation 
allows for attachment to a solid support similar to the direc- 
tion displayed as a gill fusion. Some of the peptides, such as 
R5P2 for raloxifene and 30P2 for GI5986, did not display high 
signal-to-background (at 10 nM ERa LBD) as synthetic pep- 
tides relative to signal achieved with phage ELISA. It is pos- 
sible that these peptides require some aspect of the gill 
protein or some other phage coat protein as structural scaf- 
folding to allow display of active peptide on phage, whereas 
the pure synthetic peptide is not active. Therefore, these 
weaker affinity synthetic peptides were not typically useful in 
the multiplexing format used in this report (see below). An 
alternative for future consideration is that the peptide binding 
site density on the surface of the microsphere can be tailored 
according to the affinity so that both high and lower affinity 
peptides can be used in the same multiplexed experiment 
(25). 

Fluorescent Microspheres 

Lumavidin-coated polystyrene microsphere populations (5.6 
Jim in diameter) with 10 6 biotin-binding sites per bead were 
purchased from the Luminex Corp. (Austin, TX). 

Coupling of Biotinylated Peptides to Microspheres 

Each biotinylated peptide was coupled to a specific popula- 
tion of microsphere as described previously (25). Briefly, in 
96-well 1.2-/im filter-bottom plates (Millipore, Bedford, MA), 
Lumavidin-coated microspheres were washed in buffer (PBS 



containing 0.02% Tween 20, 0.1 % BSA, 0.02% sodium azide 
and 1 mM dithiothreitol) and incubated with biotinylated pep- 
tide (2.0 /xl of a 1 00 fiM stock) for 30 min at room temperature 
in the dark. An additional 30-min incubation with free D-biotin 
(20 fil of 5 mM D-biotin) was conducted to stop the coupling 
reaction and block any unoccupied biotin-binding sites. The 
microspheres were washed twice, combined and resus- 
pended in buffer. Quantitation of the coupling reaction was 
determined by flow cytometry as described previously (25). 
All peptide couplings were evaluated to assure that peptide 
density was at least 250,000 peptide molecules per micro- 
sphere. We have determined previously (25) that this peptide 
density minimizes the variation due to avidity or rebinding 
effects. The microsphere populations were coupled on the 
day of the binding experiment. 

Coupling of ERa LBD to Fiuorochrome 

The coupling of ERa LBD to the reporter fiuorochrome Alexa 
532 (Molecular Probes, Eugene, OR) was performed essen- 
tially as described previously (25). Briefly, purified, biotinyl- 
ated recombinant ERa LBD was incubated with Strepavidin 
Alexa 532 at a 1 0:1 molar ratio of receptor to fiuorochrome for 
2.0 min. The coupling reaction was stopped by adding 20 /xl 
free D-biotin (5 mM) and incubating for 30 min at room tem- 
perature in the dark. 

Multiplexed ERa LBD Binding Assay 

All components (microspheres, ligand, Alexa 532-labeled 
ERa LBD) were combined in the wells of a 96-well microtiter 
plate in a total volume of 75 jllI. Each microtiter well contained 
3,000 microspheres of a given population and up to 60 indi- 
vidual microsphere populations. The suspensions incubated 
for 1 .5-2 h at room temperature in the dark. Fluorescence 
associated with each microsphere was directly measured by 
flow cytometry without washing. 

Microsphere-associated fluorescence was measured us- 
ing a LX-100 flow cytometer (Luminex Corp., Austin, TX) 
equipped with an XY platform for automated sampling from 
the wells of a 96-well microtiter plate. The LX-100 has two 
lasers. The red laser (633 nm) excites the fluorochromes 
embedded within the microsphere for population identifica- 
tion and the green laser (532 nm) excites the reporter fiuoro- 
chrome on the nuclear receptor to indicate binding. Instru- 
ment calibration was conducted using LX-100 calibration 
microspheres purchased from the Luminex Corp. according 
to the manufacturer's instructions. Microsphere fluorescence 
was determined per well from a minimum of 100 micro- 
spheres of each microsphere population and was gated to 
exclude doublets and aggregates using side scatter mea- 
surements. Microsphere populations were identified by their 
unique two-dimensional fluorescent profiles at 658 and 712 
nm. ERa LBD binding was measured at 575 nm and is rep- 
resented in relative MFI units. 

ERE Assay 

The ERE assay was conducted by transiently transfecting a 
T47D breast cell line with a plasmid containing 1) the frog 
vitellogenin promoter which drives expression of a secreted 
placental alkaline phosphatase reporter gene, and 2) a plasmid 
that expresses the full-length ERa gene. The vitellogenin pro- 
moter contains the consensus vitellogenin ERE sequence 
(F-vitERE) (36, 64). In this ERE assay, the presence of an ERa 
agonist increases transcriptional activity of the reporter gene. 
Compounds were evaluated for their potential to either stimulate 
transcriptional activity (agonist mode, increasing concentrations 
of compound are incubated with cells in the absence of a 
competing compound) or inhibit transcriptional activity (antag- 
onist mode, increasing concentrations of compound are incu- 
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bated with cells in the presence of a fixed concentration of 
estradiol). Transcriptional activity was measured by acquisition 
of 0D 405 which indicates secreted placental alkaline phospha- 
tase activity. An EC 50 (agonist mode) or IC 50 (antagonist mode) 
was determined for each compound tested. A compound's ERE 
potency is defined as pEC 50 x % max (agonist mode) or pICso x 
% mIn (antagonist mode). 

MCF-7 Cell Proliferation Assay 

Cell proliferation of MCF-7 breast cancer cells was deter- 
mined by [ 3 H]thymidine incorporation (37, 38). Compounds 
were tested for their potential to either stimulate MCF-7 pro- 
liferation or inhibit proliferation. An EC 50 or IC 50 was deter- 
mined for each compound tested. A compound's MCF-7 
potency is defined as pEC 50 x %cpm max (relative to estra- 
diol) or plC 50 x %cpm rT , in (relative to raloxifene) (cpm = 
radioactive counts per minute). 

Ishikawa Endometrial Cell Stimulation Assay 

The stimulation of human endometrial cells (Ishikawa cell line) 
was determined by increase in intrinsic alkaline phosphatase 
activity (39, 40). Compounds were tested for their potential to 
stimulate endometrial cells, as compared with estradiol. The 
%E max is defined as 100 x (spectrophotometric reading at 
405 nm with compound minus blank)/(spectrophotometric 
reading at 405 nm estradiol minus blank). 

Data Reduction Methods 

The multiplexed ERa LBD binding assay yielded relative flu- 
orescence intensities for 405 compounds and 58 peptides. 
Figure 2B shows these intensities (vertical axis) plotted for 
each peptide tested (horizontal axis). The points representing 
the intensities for any given compound were connected with 
line segments to yield a profile (Spotfire DecisionSite, Spot- 
fire, Inc., Somerviile, MA). This 58-dimensional data set was 
simplified by applying PCA to yield the peptides (dimensions) 
deemed to provide most of the covariance in the data set. 
PCA was performed using the JMP application (SAS Institute, 
Cary, NC), yielding eigenvectors representing the relative 
contributions of the principal components to the total covari- 
ance. The 10 eigenvectors with greatest magnitude, repre- 
senting approximately 99% of the covariance were selected. 
From these, peptides with coefficients greater than 0.25 were 
chosen for profile analysis. This resulted in a reduced data 
set, yielding profiles having 10 peptide dimensions (Fig. 2C). 

These peptide binding profiles were subjected to hierar- 
chical clustering according to profile shape. For each profile 
of N responses, all N(N-1)/2 pairwise differences were calcu- 
lated and placed into a list. The distance between two such 
lists could then be characterized by their Euclidian distance. 
Using such a metric, agglomerative clustering was applied 
(unweighted pairwise group mean average, Fig. 3). Drawing a 
boundary through the clustering tree at a level yielding 12 
clusters, and plotting profiles in a trellis arrangement by clus- 
ter ID gave a useful classification of profiles (Fig. 4). Cell assay 
data were similarly plotted as profiles. These data were lim- 
ited in availability, and so were plotted without line segments 
connecting the data points. By arranging cell assay data 
profiles in a trellis using the 12 cluster IDs, the correspon- 
dence of peptide profile shape to cell assay results can be 
visually analyzed (Fig. 8). 
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